Detecting gravitational waves from inspiraling binaries with a network of 
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We compare two strategies of multi-detector detection of compact binary inspiral signals, namely, 
the coincidence and the coherent. For simplicity we consider here two identical detectors having the 
same power spectral density of noise, that of initial LIGO, located in the same place and having 
the same orientation. We consider the cases of independent noise as well as that of correlated 
noise. The coincident strategy involves separately making two candidate event lists, one for each 
detector, and from these choosing those pairs of events from the two lists which lie within a suitable 
parameter window, which then are called as coincidence detections. The coherent strategy on the 
other hand involves combining the data phase coherently, so as to obtain a single network statistic 
which is then compared with a single threshold. Here we attempt to shed light on the question 
as to which strategy is better. We compare the performances of the two methods by plotting the 
receiver operating characteristics (ROC) for the two strategies. Several of the results are obtained 
analytically in order to gain insight. Further we perform numerical simulations in order to determine 
certain parameters in the analytic formulae and thus obtain the final complete results. We consider 
here several cases from the relatively simple to the astrophysically more relevant in order to establish 
our results. The bottom line is that the coherent strategy although more computationally expensive 
in general than the coincidence strategy, is superior to the coincidence strategy - considerably less 
false dismissal probability for the same false alarm probability in the viable false alarm regime. 

PACS numbers: 95.85.Sz,04.80.Nn,07.05.Kf,95.55.Ym 



I. INTRODUCTION 

The existence of gravitational waves (GW), pre- 
dicted in the theory of general relativity, has long 
been verified 'indirectly' through the observations 
of Hulse and Taylor. The inspiral of the mem- 
bers of the binary pulsar system named after them 
has been successfully accounted for in terms of the 
back-reaction due to the radiated GW [1]. How- 
ever, detecting such waves with man-made 'anten- 
nas' has not been possible so far. Nevertheless, 
this problem has received a lot of attention in the 
past several years, especially, due to the arrival of 
laser-interferometric detectors, which are expected 
to have sensitivities close to that required for de- 
tecting such waves [2]. The space mission LISA [3] 
is also planned by the NASA and ESA to detect 
low frequency GW. The significance of the direct 
detection of GW lies, not only in the opening of an 
entirely new window into observational astronomy; 
it further promises to place our present theories of 
gravitation on an entirely different foundation from 
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that which currently exists. Strong field gravity has 
never been directly observed and experimental ver- 
ification of the predictions of general relativity are 
very much needed. GW detection will undoubtedly 
be the most exciting development in experimental 
general relativity since the theory's birth in 1915. 

Because several detectors will take data simul- 
taneously, it is advantageous to perform a multi- 
detector search for GW signals. A multi-detector 
search would (i) improve our confidence in detection 
of a GW event if a candidate event is registered; (ii) 
provide useful directional information on the GW 
event if the detectors are sufficiently separated in 
location, that is, at large geographical separations; 
(iii) provide polarization information if they are dif- 
ferently oriented. The information mentioned in (ii) 
and (iii) is degenerate in a single detector. 

The response of a network to a GW is phase co- 
herent. By 'correcting' for the phase from the data 
of different detectors, and the time-delay between 
detectors (this amounts to bringing the detectors to 
the same site and with the same orientation), one 
can combine the data phase coherently. That this 
was achievable was shown in several works [4] . Here 
the GW source was the well studied inspiraling com- 
pact binary. 

In this paper we compare two different multi- 
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detector detection strategies the coherent versus the 
coincident. The GW source we choose is the inspi- 
raling binary because the phase of the wave can be 
sufficiently accurately computed; the phase of the 
wave has been modeled to the 3.5 post-Newtonian 
(PN) order [5] which means that the phase of the 
wave is accurate to about a cycle in a wave train 
of ~ 10 4 cycles for typical stellar mass binary in- 
spirals. Secondly, in the context of GW detectors, 
it is of great astrophysical significance [6]. Finn [7] 
has performed similar analysis with a sinusoidal sig- 
nal. However, our analysis goes further in that we 
consider a GW signal of astrophysical importance 
and our analysis uses template banks of inspirals to 
compare the two strategies. 

The coherent search strategy uses the maximum 
likelihood method where a single likelihood ratio for 
the entire network is constructed - that is the net- 
work is treated as a single detector - this is similar to 
aperture synthesis carried out, for example, by radio 
astronomers. The likelihood method combines the 
data from a network of detectors in a phase coherent 
manner to yield a single statistic which is optimal 
in the maximum likelihood sense. This statistic is 
then compared with the threshold determined by the 
false alarm rate that we are prepared to tolerate. If 
the statistic crosses the threshold, then a detection 
is announced. Note that a single likelihood ratio is 
computed and a single threshold applied in this type 
of search. 

On the other hand, the coincidence approach in- 
volves separately filtering the signal in each detec- 
tor, applying two separate thresholds corresponding 
to each detector and preparing two event lists de- 
termined by the crossings. Then the event lists are 
matched. If the estimated parameters for the events 
lie in a reasonable neighborhood in the parameter 
space of signals, a coincident detection is registered 
- if the differences in the estimated parameters lie 
within certain bounds - that is, the parameters of the 
events must lie within a certain 'window' of the pa- 
rameter space. A coincidence search for inspiraling 
binaries with real data has already been performed 
in many cases [8-10]. 

The debate between the coherent and coincident 
strategies has existed in the community ever since 
the coherent strategy was formulated. While the co- 
herent strategy involves more computational costs 
than the coincident strategy, we show in this pa- 
per for the simple case of two identical detectors in 
the same location and with same orientation that 
the performance of the coherent strategy is superior 
to that of the coincidence strategy - the false dis- 
missal is considerably less for the coherent as com- 
pared with the coincident at the same false alarm. 
This is in spite of the fact that for simplicity we 
have taken the detectors to be co-aligned; we expect 
the coherent strategy to perform better in the case 
of non-aligned detectors than its competitor. In or- 
der to decide on the performance between the two 



strategies, we plot the false dismissal versus false 
alarm curves - the Receiver Operating Characteris- 
tic (ROC) curves. We use the initial LIGO (LIGO 
I) noise curve [11] in our calculations and assume 
Gaussian stationary noise in our simulations. 

The paper is organized as follows: In section II 
we briefly review the signal in a network of detec- 
tors, the post-Newtonian inspiral waveform and also 
the maximum likelihood analysis. In section III we 
obtain analytical results of false alarm and false dis- 
missal probabilities for both coherent and coincident 
strategies. In this section we also consider the case 
of correlated noise between the detectors. This case 
would be of particular significance to the two LIGO 
detectors at Hanford and other similar future detec- 
tor topologies elsewhere envisaged. In section IV we 
perform simulations. The simulations arc performed 
in order (i) to determine the number of independent 
templates and, (ii) to determine the size of the pa- 
rameter window for the coincidence case. Simula- 
tions are necessary because at low SNR ~ 10, the 
Fisher information matrix gives poor estimates of 
the errors in parameters and therefore grossly un- 
dersizes the window. Further, simulations are also 
performed in order to plot the ROC curves which 
validate our analytical results. In section V we sum- 
marize the results and describe future directions. 



II. NETWORK SIGNAL 

A. The signal 

In this section, we describe the necessary back- 
ground and notation in order to make the paper self 
contained. Generic time domain functions h(t) will 
be denoted by h(f) in the Fourier domain, where, 

/oo 
alt h(t)e 2 * ift . (2.1) 
-oo 

Although engineering and much of LIGO analysis 
software uses the opposite sign for the exponent in 
(2.1), we maintain this notation for consistency with 
the published literature [12]. 

In order to describe the response of the detectors 
to the wave we need to define the relevant quantities 
that appear in the response. Because a lot of our re- 
sults are obtained numerically from simulations we 
give a somewhat detailed description here. The de- 
pendence of the response on the angles is best given 
by defining two frames: 

Wave frame (X, Y, Z) : The gravitational wave 
travels along the positive Z direction and X and 
Y denote the axes of polarization in such a way that 
a right-handed coordinate system is formed. 

Detector frame (x, y, z) : This denotes the orthog- 
onal coordinate frame attached to the detector. The 
arms of the detectors lie in the x — y plane, which 
is the plane tangent to the surface of the Earth and 
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the x axis bisecting the angle between the two arms 
and y is chosen such that the frame forms a right- 
handed coordinate system with the z axis pointing 
radially out of the surface of the Earth. We just 
need one such frame because the detectors have the 
same orientation. 

The Euler angles ((f), 6, ip) rotate the detector 
frame to the wave frame where we have chosen the 
Goldstein [13] convention. Also let i denote the an- 
gle between the line of sight and the orbital angular 
momentum vector of the binary. We directly write 
down the signal response h(f) in the Fourier domain. 
In the stationary phase approximation the spin-less, 
restricted post-Newtonian inspiral signal is given by, 



The function t c , 5 C , To, T3) describes the phase 
evolution of the inspiral waveform. We adopt the 
3PN formula given by 

7T 

*(/; t c , 8 C , r , t 3 ) = 2irft c - S c - - 

+ ii^ M/r 5/3 £ a ^ M f) m . ( 2 - 7 ) 



h(f) = AAx£(^,V; t )r 7/6 

x exp i$(f;t c ,S c ,T ,T S ) 



(2.2) 



where t c and 6 C are respectively the coalescence time 
and the coalescence phase of the binary. E(<j>, 8, if>; t) 
is the extended antenna pattern function [4] which 
depends on the orientation angles - 4 angles in this 
case - and is given by: 



1 



E= ) F^.VO + cos^OM.VO 

(2.3) 

where the F +iX are the usual antenna pattern func- 
tions: 



1 + cos 2 9 

F + = — sin2<^> cos2?/> — cos 20 cos # sin 2^>, 

F x = — cos 24> cos 8 cos 2tp + sin 2(f> — - - sin 2ip. 

(2.4) 

The factor Af depends on mi,m 2 the individual 
masses of the binary and the distance r to the bina- 
ries in the following way: 

where we have resorted to geometrized units of c = 
G = 1 for convenience. M = mi +m 2 is the to- 
tal mass of the binary system and rj = m\m,ij M 2 is 
the ratio of reduced mass to the total mass M. To 
ease the computation, a new set of time parameters 
{t ,t 3 } that are functions of the masses is chosen 
such that the template spacing is approximately uni- 
form over the parameter space in these parameters. 
These parameters are defined as: 
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(2.6) 



where f a is a fiducial frequency usually chosen to be 
the seismic cut-off frequency In this paper we will 
choose this frequency to be 40 Hz. 



fc=0 



where the coefficients arc as follows: 



a = 1, 
ai = 0, 
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6848 , 



4(?rM/) 



1/3' 



(2.8) 



A 



The quantity 7 is the Euler constant, 
-1987/3080 and 6 = -11831/9240. 

Although the highest PN order obtained to date 
for a two point mass inspiral is 3.5 PN order, for a 
stellar mass binary system in the typical frequency 
bandwidth of ground-based detectors from few tens 
of Hz to kHz, the 3 PN waveform is accurate to 
within a cycle in about a total of 10 4 cycles. 



B. Matched filtering and normalized templates 

It is important to distinguish between the nature 
of parameters that appear in the above form of the 
GW chirp. The two mass parameters jl = {M, r/} or 
cquivalcntly tq , t 3 are the dynamical parameters and 
determine the shape of the chirp. Exclusive to this 
set are the 'offset' parameters A = {t c , S c } which de- 
termine the duration of the chirp or the end-points 
in the time series. The latter can be quickly es- 
timated in the matched-filtering paradigm without 
having to construct template banks over them. For 
example, spectral correlators using FFTs allows us 
to estimate correlations at all time lags and thus es- 
timate t c - Similarly, using the quadrature formalism 
to construct template banks, the coalescence phase 
can be estimated analytically. However, the dynam- 
ical parameters need to be searched over by a grid of 
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templates spanning the parameter space; thus this 
is the effective parameter space. The template bank 
is set up over p = {to,T3}. The form of (2.2) allows 
us to write the explicit quadrature representation of 
the i-th template as, 

h(f;Pi,tc,8 c ) = A(s (f; pi, t c ) cos 5 C 

+ s w/2 (f;fti,t c )sm5 c ) , (2.9) 

where, 

so(/;aMc) = is w / 2 (f;Pi,t c ) • (2.10) 

The quantity A in (2.9) is the amplitude of the wave- 
form ft. This is conveniently seen by demanding that 
both the templates so and s^/2 have have unit norm; 
i.e. the scalar products (s ,s ) = (s^/2, s^/2) = 1- 
The scalar product (a, 6) of two real functions a(t) 
and b(t) is defined as 

M.J^JtfM^to ( , u) 

J Si Sh< yJ> 

where, we use the Hermitian property of Fourier 
transforms of real functions. Sh(f) is the one sided 
power spectral density (PSD) of the noise of initial 
LIGO detectors [11]. Further, fi < f < f u is taken 
to be the effective spectral window for computing 
the scalar product. The lower frequency cut-off /; 
is dependent on the sensitivity of the detector to 
seismic vibrations and is taken to be 40 Hz = f a , 
whereas the upper frequency cut-off f u is equal to 
half the full sampling frequency taken to be 2048 Hz. 
In practice, the upper cut-off frequency for normal- 
izing the templates, we choose to be 1 kHz, that is, 
we truncate the waveform at 1 kHz. The truncation 
frequency should be chosen to be the minimum of 
fi so and 1 kHz, but for the masses that we choose in 
our simulations, 7711,7712 ~ 1.4M Q , the fi so ~ 1600 
Hz which is larger than 1 kHz. 

From the definition of the scalar product Eq. 
(2.11) we deduce from Eq. (2.10) that so and 
are orthogonal, i.e. (so,s x /2) = 0. Also from Eq. 
(2.9) we see that (ft, ft) = A 2 . 

Using the normalized templates a template bank 
is generated which covers the parameter space with a 
given minimal match usually chosen to be about 0.97 
(See [14]); the correlation of a signal with a template 
must not fall below the specified minimal match frac- 
tion for at least one template in the template bank. 
The coarsest bank satisfying these conditions is con- 
structed. Having thus set up a template bank, the 
statistic p for an output signal h(t) of the interfer- 
ometer is the maximum of the signal to noise ratio 
(SNR) over all the templates (labeled by i) and the 
time-lags labeled by t c . Thus, the statistic p is given 
[30] by, 

p = max [(s , ft) 2 + 0v2,ft) 2 ] k , (2-12) 

i,t c 

which is then compared with a pre-determined 
threshold. Any crossing above this threshold is 



recorded as a candidate event. A detailed descrip- 
tion of this procedure can be found in [15-19]. 

C. Complex templates and correlations 

In view of the fact that we require just two tem- 
plates for 5 C = 0, 7r/2, namely, s and s n / 2 , it is 
found useful to combine the two templates into a 
single complex normalized template S. We then ob- 
tain complex correlations by taking scalar product 
with the data from each detector. In particular, the 
coherent statistic is constituted from these complex 
correlations [4]. In the Fourier domain the complex 
template is given by: 

S(f) = -r 7/6 exp [ t c , S c = 0, r , r 3 )] . 

(2.13) 

The normalization factor g is chosen so that (S, S) — 
2; it is given by, 

In terms of g, the amplitude of the waveform ft is 
A = AfEg/V2. 

For the two detector network, the data consists 
of two data trains, {x 1 {t)\I = 1,2; and t e [0,T]} 
where data is taken in the time interval [0, T]. As- 
suming additive noise n 1 in each detector we have: 

a; J = ft + n J , 1=1,2 (2.15) 

Note that the signal is the same in each detector 
because we have assumed identical detectors in the 
same location and with the same orientation. The 
noise is however different in each detector. But since 
we assume identical noise PSDs and also assumed 
stationarity of the noise, the noise random variables 
satisfy the statistical property: 

(n I (f)n* I (f')) = ls h (f)6(f-f), (2.16) 

where, the angular brackets denote ensemble aver- 
age. 

If the detectors have uncorrelated noise, 

(^(/KC/O) = °- 0n thc other hand > if thc 
detectors have correlated noise, then we may 
assume: 

(n 1 (f)n*2(f)) = le(f)S h (f)5(f-f), (2.17) 

where e is real and frequency dependent. In the 
time domain, we assume e as a function of \t — t'\. 
This assumption implies also stationarity in the cor- 
related noise. Thc quantity e is, in general, a func- 
tion of frequency, because different physical mecha- 
nisms generating noise are operative in different fre- 
quency regimes. Consequently, the noise correlation 
is frequency dependent. In reality, there are several 
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complications: the noise generated by a particular 
mechanism may drift in frequency with time and/or 
the e may be dependent on the amplitude of noise, 
because of non-linear effects, etc. It is not easy to 
deal with such generic cases. However, when the 
correlation of noise is not very large, it is expected 
that the simple correlation model (2.17) we adopt is 
adequate. 

If e = 0, the noise in the two detectors is uncor- 
rected; on the other hand if e = 1 then the noise is 
perfectly correlated. The cases of physical interest 
lie in between these two extremes. However, later we 
will show that we can treat e effectively as constant, 
thus simplifying the analysis. In the following, we 
analyze both cases of uncorrelated as well as corre- 
lated noise. 

We end this section by defining the complex cor- 
relations C 1 which are particularly required in the 
coherent detection strategy. We retain the notation 
and definition as in [4]. The complex conjugate of 
C 1 is denoted by Cj. It is given by: 



C; - {S,x I ) = c I -ic 



it/2 : 



(2.18) 



where, Cq and , 2 are the real and imaginary parts 

of C 1 ; they are obtained by taking the scalar prod- 
ucts of the data x 1 with s and s^/ 2 respectively; 
that is, 



C 0,tt/2 = (s(W2,2! 7 )- 



(2.19) 



III. COHERENT, COINCIDENT 
STRATEGIES AND FALSE ALARM, FALSE 
DISMISSAL RATES 

This section serves as a prelude to the next sec- 
tion, section IV, where we obtain the ROC curves 
for each strategy. As mentioned before, these curves 
are obtained by plotting the false dismissal rate (or 
equivalcntly the detection probability) versus the 
false alarm probability for a given strategy. A fair 
comparison of the strategies is obtained by plotting 
the ROC curves for each such strategy. We consider 
the following four cases: 

A. 1. Coherent detection and uncorrelated noise, 

A. 2. Coherent detection and correlated noise, 

B. 1. Coincident detection and uncorrelated noise, 
B. 2. Coincident detection and correlated noise. 

A comparison of the performance of the strategies 
is drawn pairwise between cases A. 1. and B. 1, and 
A. 2 and B. 2. 

In this section we obtain useful analytic formulae 
for the false alarm and false dismissal probabilities. 
We find that the probabilities are obtained in terms 
of certain parameters, the number of independent 
templates AT ind and the size of the parameter window 



iV w i n , which are difficult to compute analytically - we 
elaborate on these issues later in the text. However 
these quantities can be determined numerically via 
simulations. We pin down these parameters later 
in section IV so that the relevant probabilities are 
fully determined. Only then the ROC curves can be 
plotted for each strategy. 



A. Coherent detection 

In the coherent strategy, the data from differ- 
ent detectors are combined phase coherently so that 
the network operates effectively as a single detec- 
tor. The maximum likelihood network statistic for 
coincident and co-aligned detectors has already been 
worked out in [4] for the case of uncorrelated noise. 
The maximum likelihood statistics incorporates au- 
tomatically the phase of the signal and results in 
the coherent statistic. The maximum likelihood ap- 
proach tends to outperform any other adhoc test and 
we follow this approach for the case at hand. A de- 
tailed discussion on this issue, pertaining to aligned 
and misaligned detectors has been investigated by 
[20]. 



1. Uncorrelated noise 

The coherent network statistics for this case is 
given by: 



L=—\C* 1 +C* 2 \ 



(3.1) 



where the C 1 are the complex correlations of the two 
detectors / = 1,2. Although the statistic L has the 
advantage that it is proportional to the amplitude 
of the GW, the square of L has simpler form of the 
probability distribution if the noise in the detectors 
is Gaussian and stationary. Thus we define: 



A : 




(3.2) 



Equivalently, A can also be used to decide the pres- 
ence or absence of a signal. Our goal is to obtain the 
false alarm and detection probabilities Pfa and Pr> e 
respectively for the detection strategy. The false 
alarm probability is obtained in two steps: we obtain 
the probability in the single template case (for in- 
stance, when the signal parameters are known) and 
then deduce the probability where one must search 
over the parameter space. 

Consider first the single template case: 
Note that when the noise is uncorrelated, Cg and 
</y 2 , I = 1,2 are mutually independent Gaussian 
variables distributed standard normal; that each 
have mean zero and standard deviation unity. This 
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is when there is no signal present in the data. From 
this it follows that the quantities -^j(cj + Cq) and 

~^2^ c n/2 + c tt/2) are a ^ S0 tw0 dependent Gaussian 
random variables with zero mean and unit variance. 
Then A which is the sum of squares of these two 
random variables has the probability density func- 
tion (pdf) given by: 



Po(A) 



1 



-cxp 



(3.3) 



Given a threshold A*, the statistic A exceeds the 
threshold A* with the probability: 



^1 template 
FA 



/>oc 
J A' 



dAp (A) — exp 



A* 

T 



(3.4) 



Pi 



1 template 
FA 



is the false alarm probability for one tem- 
plate. Normally the threshold A* is chosen suffi- 
ciently high so that few false alarms occur during 
the observation period. 

However, the parameters of the signal are not 
known a priori and one must scan the data with a 
bank of templates. Filtering the data through all the 
templates in the bank enormously increases the false 
alarm rate. Typically, the templates densely cover 
the parameter space in order to minimize the chance 
of a signal being missed out. We have already briefly 
described such a template bank in section II. B. For 
the spinless inspiral the template bank comprises of 
templates in the parameter space spanned by the 
parameters {to, T3}. However, when determining the 
false alarm rate, the parameter t c is as important, al- 
though one uses the Fast Fourier Transform (FFT) 
algorithm to efficiently scan over the parameter t c 
and excludes it from the standard template bank; 
one effectively also has templates in t c contributing 
to the false alarm rate. If the number of templates in 
the three parameter extended bank (which includes 
t c also as a parameter besides tq, T3) are -/Vbank and if 
they all produced independent filtered outputs and 
if the false alarm probability for a single template 
is small, then the false alarm rate for the extended 
bank is just that of the single template multiplied by 
the factor -/Vbank- However, one finds that, because 
the templates are closely packed, the filtered outputs 
are strongly correlated. The dense coverage of tem- 
plates results both from the high sampling rate in t c 
as well as the high minimal match chosen. Because 
of the correlations the false alarms are substantially 
reduced. It has been shown [21] that the effect of 
the correlations is to effectively reduce the number 
of independent random variables. We find in sec- 
tion IV from simulations that this effective number 
of independent random variables -/Vj n d is much less 
than -/Vbank. Thus the false alarm probability for the 
template bank takes the form: 



FA 



iVi„dexp 



A" 



(3.5) 



When a signal h with amplitude A is present in 
the data x 1 , the means of the correlations Cq,c^ 2 
are no more zero, although they are still Gaussian 
distributed with variance unity. Now the pdf of A 
given by Eq. (3.2) is off-centered and is given by 
[4, 22]: 



Pi( A ) = 7; exp 



A + (V2A) 2 



I [(y/2A)VX] , 



(3.6) 

where, Io is the modified Bessel function of the first 
kind of order zero. For large amplitudes A » 1, 
transforming the distribution in A to one in L, one 
obtains approximately Gaussian distribution with 
mean equal to y/2 times A. This is the familiar vN 
factor of enhancement one obtains in the signal-to- 
noise ratio (SNR) with N independent detectors. (In 
fact the relevant factor 2 in Eqs.(3.2, 3.6) is changed 
to N in the case of a network of N independent co- 
aligned and co- located detectors). Since we have as- 
sumed a network of two detectors with uncorrelated 
noise this factor is just \/~2. 

By integrating pi (A) in Eq. (3.6) we can easily ob- 
tain the false dismissal or the detection probability. 
Given a threshold A* , the false dismissal probability 
is given by 



FD 



-r 

Jo 



dA Pl (A). 



(3.7) 



The detection probability Pde is just 1 — Pfd- 



2. Correlated noise 

We now consider the case when the noise in the 
two detectors is correlated. For instance, the LIGO 
detectors at Hanford are in the same location and 
also co-aligned and therefore may share some of 
the noise sources; the same physical mechanisms 
would introduce noise in the two detectors giving 
rise to correlated noise. The future planned detec- 
tor, LCGT [23], is designed such that two interfer- 
ometers with the same length are installed in the 
same vacuum system. We thus expect some amount 
of correlation in noise. Here we consider the two 
detectors which are co-aligned and co-located with 
identical PSDs and Gaussian stationary noise. We 
assume that the noise in the two detectors satisfies 
the properties given by Eqs. (2.16, 2.17). Then the 
quantities Cq v , 2 , I = 1,2 are also Gaussian random 
variables. 

We define a complex vector whose components are 
the complex correlations C 1 : 



C T = {C\C 2 ). 



(3.8) 



The superscript 'T' denotes the transpose of a ma- 
trix. The covariance matrix of Cj is given by 



TZ = (CC f ) 



2 2e 
2e 2 



(3.9) 
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where & = (C{, C|). The quantity eo is a weighted 
'average' of the noise correlation e(/) and is given 
by: 

„_4/'><Q!«fle. ,3.10) 



Sh(f) 



The off-diagonal terms in 1Z imply correlated data 
in the two detectors. Statistically independent data 
streams are obtained by forming linear combinations 
of the data streams from each detector. This proce- 
dure amounts to diagonalizing TZ. Analogous meth- 
ods were followed in [7]. We therefore obtain the 
eigenvalues and eigenvectors of TZ. The eigenvalues 
are 2(1 ± eo). We choose two eigenvectors defined by 

(V+) T = (1,1)/V2, 

(V") T = (1,-1)/V2, (3.11) 

with eigenvalues 2(1 + eo) and 2(1 — eo) respec- 
tively. We then define a orthogonal matrix U by 
U = (V + V~) which diagonalizes TZ. We define a 
vector C = UC = (C+ Cr) T where 



V2 



(C 1 ± C 2 ) . 



(3.12) 



The new complex correlations, C^, represent the so 
called pseudo-detectors. The off-diagonal terms in 
their covariance matrix are zero. Thus, we can treat 
them as independent variables. We also observe that 
the mean of C~ is zero even in the presence of a sig- 
nal. This suggests that we need to consider only 
C + to build our detection statistic. Note that C+ 
is the output of a matched filter and thus is opti- 
mal. Using the covariance matrix 1Z, it is easy to 
verify that the real and imaginary parts of C + arc 
Gaussian distributed with variance 1 + e . In the 
absence of a signal (C 1 ) = (C 2 ) = which implies 
(C + ) = 0. Defining the coherent detection statistic 
as A = |C + | 2 we obtain the distributions for A in 
the absence and presence of a signal. In absence of 
a signal we have: 



Po(A) = 



2(l + e )' CX H 2(1 + e ) 
while in the presence of a signal with amplitude A: 



A 



(3.13) 



Pi (A) 



1 



2(1 + e ) 



cxp 



(A + 2.4 2 ) 
2(1 + e ) 



x J [A 



(1 + eo) 



(3.14) 



From these probability densities we find the false 
alarm and detection probabilities. The false alarm 
probability when the data is passed through A; nc i 
independent templates is given by: 



/>oo 

P FA (A*) = N ind / dA Po (A) 



= A ind exp 



A* 



2(1 + e ) 



, (3.15) 



and the detection probability is given by: 



Pde(A*) 



dA 



A- 2(1 + e ) 
2A 



exp 



(A + 2^ 2 



2(1 + e ) 



x I [A 



(i + e )r 



(3.16) 



As expected, when eo = 1, i.e., the detectors are 
completely correlated, the expressions for false alarm 
and detection probability are consistent with those 
of a single detector. As e decreases, the false alarm 
rate falls off more quickly and the detection probabil- 
ity increases improving the performance of the net- 
work. In the other extreme limit, when eo = 0, i.e., 
the detector noises are independent, the formulae 
reduce to those given by equations (3.5) and (3.6). 



B. Coincidence detection 

Coincidence detection is a simpler strategy where 
one detector does not 'know' about the others - the 
detectors are treated as if in isolation. The data 
from the different detectors is processed separately; 
separate thresholds are set and lists of crossings of 
the threshold or candidate events are recorded as 
separate lists. The lists are then compared for con- 
sistency in the parameters of the signal. Since a 
unique GW source is supposed to be responsible for 
the signal [24], the estimated parameters of the sig- 
nal in each detector must satisfy consistency require- 
ments. Ideally, the times of arrival of the signal in 
the detectors should not exceed the light travel time 
between the detectors; in case of the spinless inspi- 
ral we consider here, in addition to the arrival times 
matching, the estimated mass parameters {to,t 3 } 
must be identical. However, because of the presence 
of noise in the detectors, the estimates of the param- 
eters may differ from their true values and lie in a 
neighbourhood of the true values. Thus when decid- 
ing the consistency of the parameters of the signal 
one must account for the difference in estimates of 
the parameters in the different detectors. Specifi- 
cally, in the case of two detectors, a detection will 
be inferred if corresponding to a crossing in the first 
detector, there is a candidate event in the second 
detector whose parameters lie in the neighbourhood 
of the parameters of the candidate event in the first 
detector. We call this neighbourhood which is a sub- 
set of the parameter space, a parameter window de- 
noted by W. For the specific situation we consider 
of identical, co-aligned and co-located detectors de- 
tecting an inspiral signal, a detection is inferred, if 
for a candidate event in detector 1 with estimated 
parameters {t\, Tq , Tg }, there is a candidate event in 
detector 2 with parameters {i 2 , Tq, t| } which lies in 
the parameter window W, a box whose vertices are 
given by {t). ± A w t c , Tq ± A w t , t£ ± A w t 3 }. Recall 
that a candidate event in each detector is decided by 
the statistic in each detector crossing its respective 
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threshold. Note that we have taken the window to 
be symmetrical about the candidate event. This as- 
sumption is not unjustified on the basis of the sim- 
ulations which we perform in the next section. 

The question arises as to how W is determined. 
One possible way is to use the Fisher information 
matrix in order to decide W. The Fisher infor- 
mation matrix gives estimates of the errors in the 
parameters due to the noise. It also suggests that 
the errors scale inversely as the SNR. So the size of 
W depends on the SNR - the higher the SNR the 
smaller the size of W, that is, the estimated values 
tend towards the true values of the parameters for 
large SNR. However, for low SNR <~ 10 or less, it 
has been shown that the Fisher information matrix 
is a poor estimator of the errors; it underestimates 
the errors by a large margin [25] (this is shown for 
the inspiral waveform up to 1 PN order). Simula- 
tions show that the actual errors can be a factor 2 
or even 3 larger than those predicted by the Fisher 
information matrix. We therefore use simulations 
(section IV) to empirically determine the parameter 
window W. We find that the window size depends 
on the SNR and reduces with increasing SNR. Since 
the parameter space is sampled discretely with tem- 
plates, W contains a finite number N w i n independent 
templates. From the foregoing discussion 7V w i n is a 
function of the SNR and decreases with increasing 
SNR. The quantity A^ w i n enters into the false alarm 
rate as we shall see below. 



1. Uncorrelated noise 

Since each detector registers events on its own, 
the detection statistic for the first detector is just 
Ai = | C 1 1 2 (without any loss of generality, we 
may assume any one of the two detectors to be 
first). Given A^d number of independent templates 
which correspond to the three dimensional parame- 
ter space, the false alarm probability in detector 1 is 
given by: 

f°° ( A* 

p fa = N ind / dAipo(Ai) = Af ind exp ( — - 

(3.17) 

Because the detectors are taken to be identical, the 
template banks in both are the same; the template 
bank depends only on the noise PSDs which are 
taken to be identical. Therefore we do not attach 
any subscript to A^d- We set the same threshold 
A* in each detector. We then draw up lists of can- 
didate events which cross the threshold A* in each 
detector. We thus obtain two lists corresponding 
to each detector. We check for consistency between 
candidate events if a candidate event in detector 1 
has at least one candidate event in detector 2 lying 
within the parameter window W centered around 
the parameters of the candidate event in detector 1. 
The size of W is determined by A* , as seen from the 



Fisher information matrix or otherwise through sim- 
ulations. Therefore N w i n is a function of A*. Also 
by our choice of the parameters {to,^}, N w [ n de- 
pends only weakly on these parameters. We ignore 
this dependence and consider A" w i n only as a func- 
tion of A*. Suppose that the number of templates 
within the parameter window is A/ w ; njt (A*). Then, 
we have N win (A*) = a(A*)A^ win! t(A*), where the fac- 
tor a must be < a(A*) < 1. 

The probability that at least one template among 
the Awi n (A*) templates registers a false alarm is 

1 - (1 - e -^*/2) N ™ {A '\ Thus the probability of 
false alarm for the two detector network is: 



Pfa(A*) = N iad e- 



-A*/2 



1- (l-e- A */2) 



Afwi»(A*) 



This unwieldy expression for the false alarm prob- 
ability for a high threshold A* > 50, (which corre- 
sponds to an one detector SNR threshold of about 
7) can be simplified to an approximate but simple 
expression: 

-Pfa (A*) ~ iV ind 7V win (A*) e - A * 

= AWi a (A*)AW(A*)e- A * -(3-19) 

The detection probability is given as follows: The 
false dismissal in each detector for a signal of ampli- 
tude A is given by: 



FD 



1 f A ' 

= 2 J dAexp 



A + A 2 



Io(AVX). 



(3.20) 

Then the detection probability of the network is just: 



DE 



(1-^ 



FD) 



(3.21) 



2. Correlated noise 



In order to compute the false alarm and false dis- 
missal probabilities for this case, we introduce a real 
correlation vector defined by, 



C T — (Cq, C*/ 2 , Cq, C^y 2 ) • 



(3.22) 



The covariance matrix of (/ = 1,2), are given 

by, 



K = (CC T ) 



1 e \ 

1 e 

e 1 

e 1 J 



(3.23) 



The procedure of the computation of the false alarm 
and the false dismissal probabilities are as follows: 
we write down the pdf of C, then integrate over the 
angular variables (f>j (as defined in Eq. (3.27)) and 
obtain the pdfs of the quantities pj = \ C 1 \ = VA~7 
in the absence and presence of the signal, namely, 
Po(pi,P2) and pi(pi,p2) respectively. From these 
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pdfs the required false alarm and false dismissal 
probabilities are easily obtained on integration. 

Before obtaining the pdf of C, we require the de- 
terminant of 1Z denoted by det 1Z: 



detK = (1 - el) 2 



and the inverse of 1Z: 
( 



R~ 



1 



Vdctft 



1 


-eo 
V 




1 





1 






-eo 

1 



(3.24) 



(3.25) 



Under the null hypothesis Hq, when, no signal 
is present, the probability distribution of these four 
quantities will be a multivariate normal distribution 
centered at the origin, i.e., 



Po(C) 



1 



cxp 



1 



(3.26) 



We now go over to the 'polar' variables pi and the 
angular variables <pi by the relations: 



Cq = pi COS 0/, 



-tt/2 



Pi sm <pi I = 1,2. 

(3.27) 

pj is the registered SNR in the I-th detector and 
is compared with the same threshold p* in each de- 
tector. We marginalize the angular variables <pi by 
integrating over them since it is the pi which are 
compared with the threshold in this type of search. 
We transform the pdf in Eq. (3.26) into the polar 
variables (pi, cf)i) and integrate over </>/. The integral 
over <pi leads again to the modified Bessel function 
I . After including the Jacobian factor, we obtain: 



Po(Pi,P2) = 



Pi Pi 

I -el 



exp 



Pi 



p\ 



2(1 -eg) 



/ e pip 2 \ 



e 

(3.28) 

Similar arguments that led to the expression in Eq. 
(3.18) for the false alarm probability in the uncor- 
rected case need to be followed also in this case, 
except that one must be careful because of the corre- 
lated noise. The expressions therefore are more com- 
plex. We will require the Bayes theorem and con- 
ditional probabilities in order to compute the false 
alarm probability. We proceed as follows: 

We define the events A and B in the (pi,p 2 ) plane 
as follows: 

A = {{pi, P2 )\pi > p* , 0<p 2 <oo}, 

B = {( Pl ,p 2 )\0<pi <oo, P2 < p*}.(3.29) 

In our computations we require the probabilities of 
the events An B and A, which are given in terms of 
the pdf Po(pi,pi) as, 

rP* r°o 

Qi = P(AnB) = / dp 2 / d Pl p {pi,pi), 
Jo J p* 

Q 2 = P(A)= / d P2 / dp lPo ( Pl , P2 ). (3.30) 
Jo J p* 



The event A represents a false alarm occurring in 
detector 1 for one template, irrespective of detector 
2 and A n B represents the event that a false alarm 
occurs in detector 1 for one template and no false 
alarm occurs in detector 2. Then by Bayes theorem, 
the conditional probability that no false alarm oc- 
curs in detector 2 given that a false alarm occurs in 
detector 1 is : 



P{B/A) 



p(A n b) 

P{A) 



Qi 
Q 2 



(3.31) 



The probability that no false alarm occurs in W 
containing N w i n independent templates in detector 
2 when a false alarm occurs in detector 1 is just 
[P(B/A)] N ™ . The probability that at least one false 
alarm occurs in detector 2 - say event C - in the pa- 
rameter window W given that a false alarm occurs 
in detector 1 is: 



P{C/A) = 1 



Qi 
Qi 



(3.32) 



The event that a false alarm occurs in detector 1 and 
at least one false alarm in W in detector 2 is A fl C. 
The probability of this event is again given by Bayes 
theorem as: 



P{A nC) = P{A)P{C/A) = Q 2 



1- ^ 



i 



(3.33) 

Recalling that there arc -/V; n d independent templates, 
the expression for the false alarm probability for co- 
incident detection for the network becomes: 



P FA (p*) =N iad Q 2 



Qi 

Q2 



(3.34) 



The above expression goes over to Eq. (3.18) of the 
uncorrelated case when e approaches zero. Also in 
the limit of large thresholds p* >> 1, we obtain 
a simple expression for the false alarm probability 
Pfa- Define the probability Q by: 



Qo 



/ d P2 / 



dpiPo(pi, Pi) ■ 



(3.35) 



Then Q 2 = Qi + Qo and as p* — ► oo, we see that 
Qo — ► 0. Then by making a simple expansion of 
Pfa in Eq. (3.34) and keeping just the linear term 
in Qo we find that, 



Pfa{p*) AWVwinQo . 



(3.36) 



This is the simple form we use in section IV to plot 
the ROC curves for which p* > 7. 

We now compute the detection probability. Under 
the alternative hypothesis Hi , when a signal of am- 
plitude A is incident on the detectors, the mean of C 
is no more zero and we have a multivariate Gaussian 
distribution translated away from the origin. Explic- 
itly, the means of its components are given by: 



(cq) = ^4cos< 



(4/2> =^4sin( 



>o , 



(3.37) 
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where 0o is the phase of the signal. Denoting the 
mean of C by C\ we obtain from the above, 



(C T ) = Cf = „4(cos 0o, sin 0o, cos 0o, sin0o) • 



Then the distribution of C is given by, 

pi(C) = 1 



(3.38) 



(27r)Vdet7e 



x exp 



(3.39) 



The exponent in pi(C) of Eq. (3.39) can be written 
as —Q/2{1 — £q), where Q is a quadratic in pj and 
is explicitly given in terms of the polar variables as: 

Q = pi + p 2 2 + 2A 2 (l-e )- 2A(l -co) 

X [pi COs(0i - 00 ) + Pi COS(0 2 - 0o)] 

- 2e O /9i/5 2 cos(0i - 2 ), (3.40) 

Integrating over the angular variables and mul- 
tiplying by the Jacobian factor of p\p 2 we obtain 
implicitly the pdf in terms of pi,p2- 



/TT fTT 
-TT J — TT 



(C). (3.41) 



We can obtain a useful approximation when the 
amplitude of the incident wave is large, i.e. 
A » 1. We evaluate the above integral using sta- 
tionary phase approximation. Since the result is in- 
dependent of 0o, the phase of the signal, we may 
without loss of generality put it equal to zero. Then 
the integrand contributes only in the neighbourhood 
of 0i = 02 = 0. Expanding Q up to second order in 
0/ we obtain: 



Q - £ 



7=1 



A2 1 



pj- 2A(l -e ) Pl (l-f ) 



+ 2.4 2 (l-£ ) 



- 2e p 1 p 2 



1 - 



{<h ~ fa? 



(3.42) 



Since the contribution to the integral over 0/ comes 
from the region mainly near 0/ = 0, we can change 
the integration region from [— 7r,7r] to [— oo, +oo]. 
We thus obtain: 



Pl(Pl,P2) 



VpTpH 



27r(l-eoM 



fo(pi + P2) 
(1-eoM 



-1/2 



exp 



(pi - A) 2 + (p 2 - A) 2 - 2e (pi - A)(p2 - A) 



2(1 -eg) 



r 



(3.43) 



This is similar to the familiar Gaussian distribution 
centered at the amplitude in both detector dimen- 
sions. Because of the correlation, the Gaussian is 
tilted at an angle of 45° with standard deviations 
yTdTeo along the two eigendirections. 

The detection probability is then given by inte- 
grating from p* to 00 along both detector dimen- 
sions. Thus, 



DE 



= / dpi \ 

J p* Jo' 



dp 2 pi{pi,p 2 )- 



(3.44) 



independent templates in all cases. These quantities 
are difficult to determine analytically and must be 
obtained from simulations. Therefore the first goal 
of this section is to estimate these quantities and in 
doing so obtain the template density in {tq,ts} for 
the 3 PN inspiral waveform. The template density is 
required to compute the number of templates iV w i n 
in W. The final goal of this section and also the pa- 
per is to plot the ROC curves for each of these cases 
and compare them for performance metrics. 



IV. NUMERICAL RESULTS FROM 
SIMULATIONS 

In the previous section we obtained useful ex- 
pressions for the false alarm and detection proba- 
bilities in the coherent and coincidence detection in 
the uncorrelated and correlated noise cases. These 
probabilities are required in the plotting of ROC 
curves in each of the four cases. However, the false 
alarm probabilities in each case Eqs.(3.5), (3.15), 
(3.18), (3.34) depend on two undetermined quanti- 
ties namely, iV w i n the size of the parameter window 
W in the coincidence cases and iV ind , the number of 



A. The template bank for the 3 PN waveform 

At high minimal matches MM close to unity, the 
template bank is obtained conveniently from the 
metric on the intrinsic parameter space {tq^t^} - 
the quadratic approximation near the peak of the 
ambiguity function is adequate. Recent work in [5] 
deals with 3.5 PN inspiral waveform and implicitly 
computes the Fisher information matrix which is in 
fact the metric on the full parameter space which 
includes other kinematical parameters on which the 
waveform depends. The idea of a metric was intro- 
duced in [25] and its computation on the space of 
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the two masses [14] - the intrinsic parameter space 
- has been amply dealt with in the literature. Here 
we do not go into the details of its computation but 
directly state the results. 

The metric depends on the noise PSD of the de- 
tectors. Since we take identical detectors the noise 



where, x = f/fk, Ik = 150Hz, f s = 40Hz, and S = 
9 x 1(T 46 /Hz. 

The metric on the parameter space {tq, t%} is de- 
noted by g a p where the indices a, (3 refer to the co- 
ordinates A" = {t , r 3 }. The significance of the met- 
ric is that the "squared distance" between adjacent 
normalized templates separated by the coordinate 
vector AA Q is given by g a pAX a AX 13 . 

This squared distance is determined by the frac- 
tional mismatch that one is prepared to tolerate i.e.: 

g af3 AX a AX 13 = 1 - MM . (4.2) 

where MM is called the minimal match. The Fig.l 
shows the contour of the templates which are at a 
fractional match of 0.97 from the central template 
(Atq = AT3 = 0) for the 3 PN inspiral waveform. 
This is in fact the contour of the ambiguity func- 
tion at the level MM = 0.97. In general the metric 
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FIG. 1: Plot of the contour of the ambiguity function at 
MM = 0.97. The fiducial frequency which scales both 
the axes is f a —40 Hz. 

would depend on the location (to, T3) in the param- 
eter space; however, the specific parameters {to,t 3 } 
have been chosen so that the metric is almost con- 
stant over the full parameter space. That this is 
possible is a consequence of the fact that the param- 
eter space is almost fiat. 



PSD is the same. We take this to be that of initial 
LIGO as given by the analytic approximation in [26] . 
Our final results are not very sensitive to the exact 
PSD and therefore it is adequate to employ the an- 
alytic approximation which is convenient to use in 
our simulations. The expression is: 



(4.1) 

I 

We may construct a template bank with the tem- 
plates placed parallel to tq and T3 axes. For MM = 
0.97, the distance Atq and AT3 between the adja- 
cent templates in each direction tq and T3 respec- 
tively is At ~ 25ms and AT3 <~ 19ms. For the 
purposes of our simulations it is sufficient to select 
a 'small' rectangle compared with the full parame- 
ter space which usually covers a region correspond- 
ing to Mi < mi,m2 < M u , where Mi < 1M and 
M u ~ 30Mq. By doing so, we have the important 
benefit of saving enormously on the computational 
costs, without sacrificing on the results as the metric 
is nearly flat. We thus assume an uniform distribu- 
tion of templates in the selected rectangle, which 
has a corner on the point (t ,t 3 ) corresponding to 
each component mass equal to 1.4M Q . We choose 
the rectangle large enough to accommodate 625 tem- 
plates. 



B. Estimating Ni n d 

Since the data is sampled at A -1 > 2048 Hz, it 
produces time samples in the statistic with sampling 
interval A < 0.5 ms. Thus the statistic is very finely 
sampled in t c which gives rise to strong correlations 
between adjacent samples. Although this is true for 
the full parameter space the correlations are partic- 
ularly strong in the parameter t c . These correlations 
tend to give us erroneous estimates of the false alarm 
probability if we consider all the samples as inde- 
pendent random variables. We need to factor in this 
effect of correlated output to estimate the number 
of independent templates M n( j. A single detector 
analysis is sufficient. 

It is known that the filtered output passed through 
a band-pass filter is correlated. It was shown in 
[27] that the chirp filter (taken to be Newtonian 
here for simplicity) acts effectively as a band-pass 
filter with the lower limit frequency as the seismic 
cut-off which acts like a 'wall' and the upper fre- 
quency determined by the signal power fall off which 
is ~ /~ 7 / 3 for the chirp waveform. Even if the raw 
noise is white, the chirp filter output for successive t c 
is correlated, because the band-pass filter effectively 



J 



Sh(f) = S [(4.49a;)- 56 + 0.16a;- 4 - 52 + 0.52 + 0.32a; 2 ] / > f s 
= 00 otherwise, 
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reduces the number of degrees of freedom in the orig- 
inal raw data. In the figure below we plot the auto- 
correlation function of the filtered output c(t c ) as a 
function of At c , namely, a(At c ) = (c(t c )c(t c + At c )) 
which is maximized over the initial phase but with 
the same mass parameters. The At c at which the au- 
tocorrelation approaches zero is the required decor- 
relation time. From the Fig. (2) it appears that the 
decorrelation time ~ 15ms when the autocorrelation 
function falls to few percent of its maximum value. 




FIG. 2: Plots of the auto-correlation functions c(At c ) 
(dotted curve) and co(At c ) (solid curve). For reference, 
the autocorrelation curve for one of the quadrature com- 
ponents Co is also plotted. This curve is obtained both 
analytically in [27] and here numerically through simu- 
lations. The curve for c is obtained purely numerically 
via simulations. 

More directly, relevant to the false alarm prob- 
ability computation, we evaluate the decorrelation 
time of the filtered output by the following proce- 
dure. We take data trains 32 sec long which are sam- 
pled at 2048 Hz giving 65536 points in a data train. 
We take 1500 realizations of the noise and pass the 
data through 500 templates in the rectangle. We di- 
vide the matched filter output into data chunks of 
various lengths starting from the sampling interval 
of 0.488 ms up to 512 ms increasing each time the 
length of the data chunk. The size of the chunk is 
the trial decorrelation time. For each fixed length of 
the chunk (or trial decorrelation time) we maximize 
the statistic Ai = | C 1 1 2 over each chunk and the 500 
templates, and compare it with a threshold A*. We 
count the number of crossings Nfa of the statistic 
Ai > A* - the number of false alarms. We take the 
range < A* < 50 and obtain Nfa 

function 

of A*. The number of false alarms can range up to 
32 x 2048 x 1500 ~ 10 8 for the case At c = 0.488ms. 
We then plot the curve Nfa 

M.S M function of A* for 
a fixed data chunk of length At c . We now increase 
the length of the chunk At c from the sampling rate 
of ~ 0.488 ms to 512 ms successively and plot the 
curves for various chunk lengths (trial decorrelation 
times) At c . We find that for the decorrelation times 
At c > At coo ~ 15.8 ms the curves start to coalesce 



into a single curve in the large A* regime. This sig- 
nifies that the samples spaced At c > At coo apart 
are essentially decorrelated and will produce satis- 
factory results. This procedure is depicted in Fig. 3 
by plotting some of the curves with In Nfa versus 
A* for various fixed values of At c . 
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FIG. 3: The number of false alarms Npa as a function 
of the threshold A* for various fixed values of At c . 

However, in order to obtain good statistical re- 
sults, it is desirable to choose the maximum number 
of data chunks. We satisfy both these criteria by 
choosing the data chunk length to be the least, but 
at the same time ensuring decorrelation in samples, 
namely, we choose the chunk length to be At coo . 
From the simulations we have also verified that the 
autocorrelation function a(At c ) falls off to a value 
of about 3% of its maximum value at At c = 15.8 ms 
(the figure (2) also suggests this fact). 

A similar procedure is followed for evaluation of 
the decorrelation length in the parameters To and 
t 3 . We find that the samples spaced Ar > 220ms 
apart for r and At 3 > 80ms for r 3 , are decorrelated 
and will produce satisfactory results. The template 
space for our simulation has the length 600ms for 
to and 450ms for t 3 . We then take the maximum 
among the samples with different To, t 3 parameters. 

The number of false alarms at sufficiently high 
values of the threshold are given by: 



N FA (A > A*)) = N sim N ind exp 




where, N sim is the number of simulations we per- 
form. We plot In Nfa versus A* in Fig. 4 below. In 
this plot, we perform iV s i m = 1500 simulations. The 
data length is 32 seconds with a sampling rate of 
2048Hz. A rectangle with 50 x 10 = 500 templates 
was taken in the (to, t 3 ) plane. 

We observe that the graph of In Nfa versus A* is 
a straight line when A* is sufficiently large. We use 
the least square fit in order to fit a straight line in 
the regime of large A*, namely, in the range 17 < 
A* < 37. The equation of this straight line is: 

In N FA = -0.46A* + 22.96. (4.4) 
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FIG. 4: The natural log of the number of false alarms is 
plotted versus the threshold A*. In the large threshold 
regime, a straight line is fitted to this curve and extended 
to obtain intercepts on both axes. The intercept on the 
vertical axis is ~ 22.96. 



Then iV; n d can be determined by extending the 
straight line to obtain intercepts on cither axis. We 
find that the intercept on the In TVp^-axis is at 
IiiNfa = 22.96. From Eq. (4.3), we conclude that 
Wind ~ e 22 96 /1500 - 6.24 x 10 6 which is about 19% 
of the total number of templates 32 x 2048 x 500 = 
3.28 x 10 7 spanning the three dimensional parameter 
space {t c ,T ,T 3 }. 

C. Determining the parameter window W 

In this section we discuss the size of the 
parameter-window which is required to decide co- 
incident detection. Let a signal be detected individ- 
ually by the two detectors. Even though the detec- 
tors are assumed to be identical with respect to the 
noise PSD, location and orientation, the different 
noise realizations in the two detectors will in general 
produce different estimates of the true parameters 
of the signal [31]. If the SNR is sufficiently high 
~ 10 then presumably the detected parameters will 
not differ too much from the true parameters of the 
signal, and thus also not too much from each other. 
We therefore need to decide on the size of a box or 
a window in the parameter space, in order to say 
whether it is the same signal from an unique astro- 
physical source detected in each detector. Taking 
into account the reasons mentioned in section III B, 
we determine the size of W by performing simula- 
tions. 

As seen from the arguments pertaining to the 
Fisher information matrix or from the simulation 
results, the size of W depends on the SNR of the 
signal - the higher the SNR, the smaller the size of 
W - the error in the estimates reduces. For a signal 
injected with typical parameters (lying somewhere 
near the center of the parameter space) with typical 
SNR = 10, the Fig. 5 shows the distribution of the 



difference {At c , Ato, AT3} of the parameters in two 
detectors. 
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FIG. 5: Distribution of the difference of the detected 
parameters, (a) t c , (b) to and (c) T3, for an injected 
signal of SNR=10. 

750 simulations have been performed for two de- 
tectors and the results are plotted as histograms. 
The vertical axis shows the number of events de- 
tected in each bin for each of the parameters in the 
three figures. The distributions taper off as we go 
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away from the true signal parameters. We fix the 
size of W by requiring that at least 99% of the ob- 
served events in each parameter fall in the intervals 
|Ai c | < A w t c , |At | < A w t and |Ar 3 | < A w t 3 . 

I 



W = {(At c ,Aro,AT 3 )||At c | < A w t 



The above results were obtained for a typical SNR 
of 10. Because of noise, even if the signal has true 
SNR of 10, the observed SNR may differ. For a 
signal with true SNR = 10, the observed SNRs are 
found to vary anywhere from roughly 7 to 14. This 
is shown in Fig. (6) below. 
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FIG. 6: Distribution of SNR for an injected signal of 
SNR=10. 

However, W depends on the SNR. So for the three 
parameters we determine W for various SNRs rang- 
ing from 7 to 14 in steps of unity. The window- 
sizes were taken such that at least 99% of the ob- 
served values in each dimension fall within the win- 
dow. This procedure detects the signal on an average 
(0.99) 3 ~ 0.97 of the time. 

The following table (I) shows the results of our 
simulations. 750 signals of a specific SNR, say 7, are 
injected in the simulated noise of two detectors and 
then the number of events within a given bin size 
are counted and this number is normalized to unity 
by dividing by 750. This simulation obtains for us a 
window-size in each parameter (three in all for the 
three parameters) for SNR = 7. The procedure is 
repeated in addition for SNRs ranging from 8 to 14 
in steps of unity. From the table we may infer the 
window size in each parameter for the various SNRs. 
The window size is tabulated in units of the bin size, 
namely, the distance between templates, for each of 
the parameters. We err on the safer side by choosing 
the bin size to be the smallest integral value larger 
than the corresponding fractional bin size for which 



This requirement fixes the window size. The per- 
centage 99% is subjective but seems to us to be rea- 
sonable. We therefore have: 



, | Aro | < A w r , |Ar 3 | < A w r 3 } (4.5) 
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TABLE I: Window size in the parameters t c , to and T3 
for SNRs ranging from 7 to 14 in steps of unity. Window 
size corresponds to 99% detection. 



the cumulative frequency crosses the 0.99 mark. 

The window size W determines the number of 
template N w - ln ^ t within the windows. In order to 
plot the ROC curves, we need to estimate N win , or 
equivalently a. However we find that it is not easy 
to determine a as we did to determine N ind . We 
thus need the help of the full numerical simulation 
to determine a, which is discussed in the next sub- 
section. 



D. The ROC curves 

As mentioned before, the ROC curves fairly com- 
pare the coherent and coincidence strategies of de- 
tecting signals. At the same false alarm we can com- 
pare the detection efficiency (probability) by plot- 
ting these curves. We consider several situations 
here. The simplest of these involve injecting signals 
in a small rectangle in the (to, T3) space in the neigh- 
bourhood around the point corresponding to indi- 
vidual masses of 1.4M . We plot the ROC curves 
for the uncorrelated case eo = and the correlated 
cases eo = 0.2, 0.3. We then extrapolate our results 
to the full parameter space by making certain as- 
sumptions which we state later and plot the ROC 
curves. Finally, we consider a uniform distribution 
of sources, uniform in spatial density and orientation 
with individual masses distributed close to 1.4M Q 
and compare the performance of the two strategies. 

The signal is injected in a rectangle IZs of the 
( t 0j7"3) plane with its center at tq = 25.1 sec , 
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T3 = 1.0743 sec which lies well inside the parame- 
ter space, and away from its boundaries. The sig- 
nal parameters are drawn randomly from a uniform 
distribution over the section of the parameter space 
consisting of the rectangle IZs in the (rcTa) space, 
t c and initial phase. The SNR is fixed at 10. In 
order to detect a signal in 7?..$, the detection tem- 
plates also span a rectangle TZr, larger than IZs to 
take into account effects of noise - a signal lying in 
IZs could be detected by a template outside it, es- 
pecially if the signal lies close to its boundary. TZq- 
is chosen sufficiently large so that the width of the 
margin TZt ~ TZs is ~ ^-w T o an d ~ ^w T 3 m the 
relevant dimensions. 

Specifically, we take IZs = (24.9375 — * 
25.2625) sec x (0.9603 -> 1.2073) sec and 1Z T = 
(24.8 -> 25.4) sec x (0.8463 -> 1.3023) sec. 

The resulting figures are plotted for eo = 0, 0.2 
and 0.3 (figures 7, 8(a) and 8(b), respectively). The 
trend of declining performance for increasing val- 
ues of eo is evident. For the sake of comparison, 
we have also plotted the ROC curve of a single de- 
tector. However, the performance of coherent de- 
tection is better than that of the coincident detec- 
tion for all thresholds. We get a factor of roughly 
2 to 4 times improvement in the detection efficiency 
over the coincident detector for a given rate of false 
alarm. It justifies, why despite the heavier compu- 
tational burden [24] , the coherent detection strategy 
is preferable. It allows one to push the false alarm 
rate down for a given detection efficiency, increasing 
the confidence of detection. In Figs. 7 and 8, the- 
oretical ROC curves derived in Section III are also 
plotted. We use the value of A^d determined in 
Section IV. B. On the other hand, we set a = 1 (i.e., 
N w - ln = N w - m , t ) for the value of window parameters. 
We find that with this choice of the value N- m< i and a, 
we can obtain a fairly good fit of theoretical curves 
to the simulation curves. As discussed in Section 
IV. C, the size of parameter windows and the decor- 
relation length suggests that the number of indepen- 
dent templates in the parameter window is nearly 1 
(i.e., a ~ 1/-/V w i njt ). Although a small change in a 
does not change the ROC curve very much, if we 
set a ~ l/iV w i ni t, the theoretical ROC curve devi- 
ates from the simulation curve significantly. This 
suggests that although the templates within a pa- 
rameter window in Table I are spaced within the 
decorrelation length, the correlations between sam- 
ples are not very large. We thus need to take a 
which is slightly smaller than unity, a < 1. 

Although, it is not very easy to explain the reason 
for this choice of a for coincident detection, since the 
theoretical ROC curves and simulated ROC curves 
agree very well, we can use the theoretical curves for 
the evaluation of the performance of the coherent 
and coincident search of a much larger data set. We 
consider a one year data train and search the whole 
parameter space spanned by masses 1M© - 40M©. 
We assume that the sampling rate of the detector is 
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FIG. 7: The ROC curves for single detector, uncorrelated 
coincident and coherent analysis for an injected signal of 
SNR=10. The solid line, dashed line and the dotted line 
correspond to the theoretical ROC curves for coherent, 
coincident and single detectors derived in Section III, 
with iVjnd determined in Section IV. B. We also take a = 
1 irrespective of the threshold. 



2048Hz. Then the number of time templates would 
be 

length of interval 3 x 10 7 sec „ „ „ 1f 

N t = — -. ~ 6.1xl0 lc 

sampling interval 4.9 x 10 4 sec 

(4.6) 

In order to evaluate the number of mass templates, 
we assume the metric is constant in the parameter 
space, for simplicity. We follow the analysis by Owen 
[28]. For the LIGO I noise spectrum, 3PN phasing 
formula and MM = 0.97, we obtain V m ~ 1.2 x 10 4 . 
Thus, the total number of templates is 

N tot = N t N m ~ 7.3 x 10 14 . (4.7) 

In section IV. B, we have determined the value of 
this fraction of independent templates to the number 
of total templates to be 0.19. As only the parame- 
ter space has enlarged and the metric is nearly flat 
(which implies the correlations between templates 
depend essentially on the difference between param- 
eters and not on the locations of the templates), the 
fraction of independent templates to the total num- 
ber of templates should not be very different from 
this value - we assume this fraction to be the same, 
namely, 0.19. We also assume a = 1 for the co- 
incident detection. Assuming Vi nc j and a in this 
manner, and SNR = 10, we plot the ROC curves in 
Fig. 9. We find that, in fact, the net effect is to 
only translate the ROC curves along the Pfa ax is 
relative to the ones plotted in the previous cases. 

Finally, we consider the case of inspiral sources 
distributed uniformly in space and uniformly in ori- 
entation - the direction of the orbital angular mo- 
mentum of the binary is distributed uniformly over 
the sphere. We assume a uniform distribution within 
0.4 Mpc and 26 Mpc. The lower limit on the dis- 
tance ensures that our galaxy is excluded and the 
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FIG. 8: The ROC curves for coherent and coincidence 
detection for two correlated detectors with eo values (a) 
0.2 and (b) 0.3 for an injected signal of SNR=10. 
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FIG. 9: ROC curves for single (dotted line) coher- 
ent (solid line) and coincident (dashed line) detectors (un- 
corrected) for a 1 year data train sampled at 2048Hz 
and searched over IMq — 40Mq, for an injected signal 
of SNR=10. 

upper limit on the distance comes from the maxi- 
mum range of LIGO I detectors for SNR of 7. The 
direction of the orbital angular momentum is de- 
termined randomly. We consider the one year data 



train, sampled at 2048Hz and the mass parameter 
region, 1 — 4OM . The ROC curves are given by 
Fig. 10. We find that the coherent detection is bet- 
ter than the coincident detection by about 25-40%. 
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FIG. 10: ROC curves for single(dotted line) coher- 
ent (solid line) and coincident (dashed line) detectors (un- 
corrected) for a 1 year data train. The sources are dis- 
tributed uniformly in distance between 0.4 Mpc and 26 
Mpc and also in orientation. 



V. CONCLUSION 

We have compared two possible strategies of 
multi-detector detection, namely, the coincidence 
and the coherent. We have considered the inspiral 
waveform for this purpose and by plotting the ROC 
curves for the two strategies we obtain a fair as- 
sessment of their performance. We have considered 
several cases from the simple to implement, to the 
more astrophysically relevant one. The bottom line 
here is that the coherent strategy of detection is su- 
perior to the coincidence strategy. The ROC curves 
display this fact quantitatively for all the cases con- 
sidered. The coherent detection strategy uses the 
likelihood function which inherently uses the infor- 
mation of phase coherence to decide on detection. 
This information is encoded in the pdf of the H\ hy- 
pothesis, that a signal is present in the data. The 
pdf explicitly contains the signal from each detector 
in the network with consistent phase. Thus the co- 
herent strategy inherently accommodates the phase 
information. Moreover, the likelihood analysis leads 
naturally to the matched filter which is the optimal 
filter for the network output. On the other hand, the 
coincidence strategy uses separate event lists formed 
by identifying candidate events in each detector con- 
sidered in isolation - it ignores the crucial phase in- 
formation that is inherent in a signal generated from 
a specific astrophysical source. This critical differ- 
ence between the two analyses leads to the coherent 
analysis having superior detection performance. 

Moreover, when we consider correlated noise be- 
tween detectors, the coherent strategy again leads 
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to matched filters which are optimal. The correla- 
tions between the noise enter naturally into the like- 
lihood analysis and the matched filter depends on 
these correlations. In coincidence analysis, on the 
other hand, the correlations between detectors sim- 
ply do not play any role in the detection procedure 
and are thus ipso facto ignored. 

In this paper we have treated two detectors with 
identical noise PSD, identically oriented and in the 
same location. The results are easily generalized 
from two to N detectors with identical noise PSD, 
with identical orientation and in the same location. 
We next propose to generalize these results to in- 
clude arbitrarily oriented detectors in different geo- 
graphical locations. We expect that the difference in 
performance between the two strategies to be more 
striking. One reason is that for differently oriented 
detectors, there is a finite possibility that a signal 
may not be detected separately by both detectors, 
thus ruling out coincidence detection, while coher- 
ently it could still be detected. 

The work in this paper can be generalized in a 
straight forward manner to other GW sources with 
known waveforms of finite duration. More impor- 



tantly, we believe that this work could be general- 
ized to burst sources where the waveform could be 
constrained from physical and astrophysical consid- 
erations. Some work already exists in this direction 
[29]. This may be a worthwhile future direction to 
adopt. 
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